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INTRODUCTION 

The  disposal  of  the  demilitarization  stockpile — unwanted  munitions,  rocket  propellants, 
and  manufacturing  wastes — is  necessary  at  Department  of  Defense  (DOD)  and 
Department  of  Energy  (DOE)  installations.  The  disposal  methodologies  include: 

1)  recovery  and  reclamation  technology,  2)  thermal  destruction  methods  such  as 
incineration  and  popping  furnaces,  3)  research  stage  technology  such  as  electrochemical 
reduction  and  biodegradation,  and  4)  open  burning  (OB)  or  open  detonation  (OD) 

(Ref.  1).  OB/OD  takes  place  in  an  earthen  pit,  trench,  or  bermed  area  and  is  the  most 
common  disposal  method  in  use  today;  this  stems  from  its  low  cost,  effectiveness,  and 
the  capacity  to  treat  a  wide  range  of  munitions. 

The  existing  demilitarization  stockpile  is  estimated  to  be  about  400,000  tons  and  is 
increasing  at  the  rate  of  about  40,000  tons  per  year.2  However,  the  material  destroyed  in 
a  single  detonation  typically  ranges  only  from  100  to  5000  lbs,  while  the  quantity  treated 
in  a  burn  is  somewhat  larger  and  usually  lasts  from  1  to  5  min.  Thus,  a  large  number  of 
detonations  or  burns  will  be  required  to  significantly  reduce  the  existing  stockpile. 

OB/OD  operations  generate  air  pollutants  and  require  predictions  of  pollutant  concen¬ 
trations  to  assess  air  quality  impacts  and  health  risks.  The  pollutants  include  SO2,  N02, 
CO,  particulates,  volatile  organic  compounds  and  hazardous  or  toxic  materials  such  as 
metals,  cyanides,  semivolatile  organics,  etc.2’3  For  very  large  detonations  (1  -  3xl04  lbs), 
natural  dust  entrained  by  the  blast  is  an  additional  contaminant  to  consider.  Emissions 
from  OB/OD  sources  have  the  following  special  features:  1)  “instantaneous”  or  short- 
duration  releases  of  buoyant  material,  2)  considerable  variability  in  the  initial  cloud 
size,  shape,  and  height,  and  3)  ambient  exposure  times  for  individual  clouds  that  are 
significantly  less  than  the  typical  averaging  times  (  >  1  hr)  of  air  quality  standards. 

Predictions  of  air  quality  impact  require  the  use  of  atmospheric  dispersion  models 
together  with  model  inputs  on  source  and  meteorological  conditions.  Currently,  there 
is  no  recommended  EPA  model  to  handle  the  special  features  of  OB/OD  sources.  The 
most  commonly-used  approach  is  INPUFF,4’5  a  Gaussian  puff  model.  The  basic  puff 
framework  is  suitable  for  OB/OD  releases  although  the  existing  INPUFF  has  several 
limitations  as  discussed  below.  As  a  result,  a  model  development  program  was  initiated 
under  the  sponsorship  of  the  DOD/DOE  Strategic  Environmental  Research  Development 
Program. 

In  the  following,  we  discuss:  1)  background  issues  influencing  the  development  of  an 
OB/OD  dispersion  model,  2)  a  model  development  overview,  and  3)  the  framework  for 
short-range  modeling  (distances  <  30  km).  Plans  for  long-range  modeling  (distances 
>  30  km)  are  in  the  initial  stages  of  development  and  will  be  described  later.  The  model 
development  program  began  in  September  1994  along  with  a  parallel  effort  to  construct 
a  mobile  meteorological  platform,  which  is  necessary  due  to  the  remoteness  of  many 
of  the  DOD  facilities.  A  related  program  has  been  acquiring  information  on  OB/OD 
emission  factors  from  experimental  test  chambers3  and  field  studies.2 


2 


95-MP22A. 05 


BACKGROUND 

Several  factors  have  motivated  and  influenced  the  development  of  an  OB/OD  dispersion 
model  including:  1)  the  limitations  of  existing  models,  2)  the  improved  knowledge  of 
the  planetary  boundary  layer  (PBL),  3)  potential  future  OB/OD  operations,  and  4)  the 
development  of  a  mobile  meteorological  platform.  These  topics  are  briefly  discussed  in 
the  following. 

Limitations  of  Existing  Models 

As  noted  earlier,  the  INPUFF  Model4-5  is  a  commonly-used  approach  for  dealing  with 
OB/OD  sources  and  can  handle  dispersion  from  individual  puffs  or  clouds  or  from 
a  sequence  of  puffs  as  in  a  short-duration  release,  e.g.,  an  open  burn.  Although  the 
Gaussian  puff  approach  is  appropriate  for  OB/OD  sources,  INPUFF  has  the  following 
limitations: 

1)  It  uses  dispersion  parameters  (<Ty ,  <jz )  from  the  Pasquill- Gifford  (PG)  curves  , 
which  are  only  applicable  to  surface  releases,  or  from  Irwin  s  scheme. 

2)  It  adopts  Briggs’  plume  rise  expressions8  which  apply  to  continuous  releases 
rather  than  to  instantaneous  sources  (puffs,  clouds,  or  thermals)  and  does  not 
address  buoyant  thermal  penetration  of  elevated  inversions  capping  the  PBL. 
Thermal  penetration  of  the  inversion  may  be  important  for  large  detonations  or 
burns. 

3)  It  assumes  Gaussian  statistics  for  turbulent  lateral  and  vertical  velocities  in  the 
PBL,  whereas  the  vertical  velocity  statistics  in  the  unstable  PBL  are  positively 
skewed.9  The  skewness  should  be  included  for  vertical  dispersion. 

4)  It  does  not  address  transport  and  dispersion  in  the  vicinity  of  shorelines, 
mountains,  and  other  complex  terrain. 

From  a  scientific  viewpoint,  use  of  the  PG  curves  is  deficient  in  that  1)  they  are  based 
on  dispersion  from  a  ground-level  source  and  for  short  downwind  distances  (<  1  km)  and 
2)  the  curve  selection  scheme  is  based  on  surface  meteorology,  which  does  not  account 
for  the  vertical  structure  of  PBL  turbulence.19  For  large  detonations  or  burns,  source 
buoyancy  can  carry  emissions  to  several  hundred  meters  or  to  the  top  of  the  PBL,  with 
the  possible  penetration  of  the  capping  inversion.  One  must  then  deal  with  dispersion 
throughout  the  entire  PBL  and  have  a  better  characterization  of  buoyancy  effects. 

Other  dispersion  models  for  OB/OD  sources  have  been  proposed  and  are  described  in 
Ref.  11. 

Turbulence  and  Dispersion  in  the  Planetary  Boundary  Layer 
Over  the  past  two  decades,  much  progress  has  occurred  in  our  knowledge  of  turbulence 
and  dispersion  in  the  PBL,12  both  for  the  unstable  or  convective  boundary  layer 
(CBL)  and  the  stable  boundary  layer  (SBL).  For  the  CBL,  numerical  and  laboratory 
simulations  and  field  observations  revealed  the  large-scale  flow  structures  and  the 
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important  turbulence  velocity  and  length  scales — the  convective  velocity  scale  w+ 
and  the  CBL  height  h .  Typical  values  of  and  h  at  midday  are  1-2  m/s  and  1500 
m.  Major  insights  into  dispersion  followed  from  laboratory  experiments,  numerical 
simulations,  and  field  observations  for  both  nonbuoyant  and  buoyant  plumes.10’13’14 

For  the  SBL,  the  turbulence  is  much  weaker  with  typical  eddy  sizes  on  the  order  of  tens 
of  meters  or  less.9  Numerical  models  and  field  observations  have  demonstrated  that 
wind  shear  is  the  important  source  of  turbulence  with  the  friction  velocity  u+  being 
the  relevant  velocity  scale;  is  typically  of  the  order  of  0.1  m/s  in  strongly  stable 
conditions.  Dispersion  has  been  put  into  a  sound  framework  for  near-surface  sources, 
whereas  the  framework  is  less  general  for  elevated  releases.15  Nevertheless,  models  have 
led  to  a  good  understanding  and  organization  of  observations. 

The  application  of  the  improved  knowledge  of  the  PBL  has  been  discussed  in  a 
number  of  short  courses  and  monographs  and  is  now  being  incorporated  into  models 
for  applications.12  A  recent  example  is  AERMOD16  for  industrial  source  complexes. 
Knowledge  of  flows,  dispersion,  and  other  processes  over  complex  terrain  is  summarized 
in  another  recent  monograph.17 

Potential  Future  OB/OD  Operations 

In  contemplating  a  significant  reduction  of  the  demilitarizion  stockpile,  consideration 
is  being  giyen  to  much  larger  detonations  (e.g.,  1-5  xlO4  lbs)  than  those  currently 
used  because  of  the  higher  temperatures  and  more  efficient  thermal  destruction  of 
contaminants.  Possible  dispersion  scenarios  include:  1)  a  large  daytime  release  with 
sufficient  source  buoyancy  to  carry  material  to  the  top  of  the  CBL  with  possible 
penetration  of  the  elevated  inversion,  and  2)  a  large  nighttime  release  with  sufficient 
buoyancy  to  carry  the  emissions  above  the  SBL  into  the  overlying  weakly-  or  non- 
turbulent  airflow.  In  both  scenarios,  the  source  material  could  be  transported  large 
distances  (~  20  to  100  km)  with  significant  lateral  dispersion  but  minimal  vertical 
dispersion,  thus  preventing  high  ground-level  concentrations  (GLCs)  near  the  source. 

The  avoidance  of  a  high  near-source  impact  would  increase  the  importance  of  long-range 
transport  with  somewhat  lower  GLCs. 

Mobile  Meteorological  Platform 

Due  to  the  remoteness  of  many  DOD  facilitites,  a  mobile  meteorological  platform 
is  being  developed  to  provide  the  PBL  variables  necessary  for  modeling.  The  initial 
platform  design  contains:  1)  a  radar  wind  profiler  for  obtaining  profiles  of  the  three 
wind  components  to  a  height  of  3  km,  2)  a  radio  acoustic  sounding  system  (RASS)  for 
temperature  profile  measurements,  3)  a  mini-SODAR  for  measuring  profiles  of  wind  and 
the  vertical  turbulence  component  ( a to  a  height  of  ~  200  m,  4)  a  mini-lidar  system 
for  measuring  the  PBL  depth,  and  5)  a  portable  meteorological  station  for  measuring 
near-surface  winds,  temperature,  turbulence,  and  heat  flux.  The  dispersion  model  should 
be  designed  to  maximize  the  use  of  the  data  from  this  platform,  with  the  temporal  and 
spatial  resolution  of  the  measurements  being  determined  by  instrument  limitations  and 
modeling  needs. 
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MODEL  DESIGN  OVERVIEW 

In  the  following,  we  give  a  brief  overview  of  the  key  features  to  be  included  in  the 
OB/OD  dispersion  model  and  the  division  by  model  types.  The  important  features  to 
address  in  the  modeling  are: 

1)  all  aspects  of  the  source  including  the  instantaneous  nature  of  the  release,  the 
cloud  or  thermal  rise,  thermal  penetration  of  an  elevated  inversion,  and  the  short 
exposure  time  of  the  cloud, 

2)  modern  dispersion  concepts12  based  on  the  turbulence  structure  and  scaling  of 
the  CBL  and  SBL, 

3)  use  of  micrometeorological  variables  along  with  vertical  profiles  of  wind, 
temperature,  and  turbulence  from  the  mobile  meteorological  platform, 

4)  short-  and  long-range  dispersion  where  the  distinction  between  them  is  taken  at  a 
distance  of  ~  20  to  30  km,  and 

5)  a  treatment  of  complex  terrain  which  exists  in  the  vicinity  of  many  DOD  facilities 
in  the  western  US. 

In  addition,  the  design  should  consider:  1)  modeling  of  the  dose  (tune  integral  of  the 
.concentration)  as  well  as  the  concentration  with  provisions  for  determining  the  time- 
averaged  concentration  that  is  necessary  in  air  quality  assessments,  and  2)  short-term 
fluctuations  in  concentration  and  dose,  and  3)  deposition  of  particles. 

The  modeling  is  divided  in  two  ways:  1)  short-  and  long-range  dispersion,  and  2) 
modeling  methodology  which  refers  to  the  degree  of  detail,  spatial  resolution,  and 
computation.  The  division  at  a  scale  of  20  to  30  km  is  somewhat  arbitrary  but  intended 
to  distinguish  a  regime  where  simple  wind  field  modeling  may  be  accomplished  (short 
range)  and  one  where  a  more  complete  wind  field  model  is  necessary  (long  range),  i.e., 
for  transport  times  exceeding  ~  1  hr.  Ultimately,  the  short-range  model  would  be  a 
component  of  or  treated  as  an  initial  “subgrid”  approach  in  the  long-range  model. 

The  modeling  methodology  is  divided  into  an  applications  approach  with  relatively 
low  computational  costs  and  a  research  model.'  For  the  applications  methodology,  a 
Gaussian  puff  model  is  proposed  whereas  a  Lagrangian  particle  model  is  planned  for 
the  research  approach.  The  applications  model  would  be  useful  for  routine  problems 
in  regulatory  permitting,  whereas  the  research  model  is  necessary  to  address  more 
complicated  issues,  e.g.,  inversion  penetration,  complex  terrain,  and  those  associated 
with  larger  detonations. 
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SHORT-RANGE  MODEL 


Applications  Model 

The  following  pertains  to  the  model  for  an  instantaneous  release  or  detonation.  The 
concentration  field  for  an  open  burn  or  short-duration  release  is  obtained  by  integrating 
the  concentration  expression  for  the  instantaneous  source  or  puff  over  time,  i.e., 
integrating  the  concentration  over  a  sequence  of  puffs  from  successive  release  times.  This 
is  briefly  discussed  under  cloud  rise  below. 


Concentration  field.  Dispersion  models  predict  the  ensemble-mean  concentration  C 
for  a  given  set  of  source  and  meteorological  conditions,  i.e.,  the  concentration  that  would 
be  observed  if  the  same  experiment — same  source  and  meteorological  conditions — were 
repeated  a  large  number  of  times.  For  now,  our  focus  is  on  the  C  for  a  given  averaging 
time,  but  it  should  also  be  possible  to  model  the  rms  concentration  fluctuation  (e.g.,  see 
Ref.  18).  In  this  section,  we  discuss  near-instantaneous  or  short-term  concentrations; 
time-averaged  concentrations  are  considered  under  the  dosage. 

Currently,  it  is  not  clear  what  short-term  concentrations  are  relevant  for  permitting 
situations  and  we  consider  two  estimates:  the  peak  and  the  mean  at  a  downwind 
location.  Both  cure  obtained  from  a  Gaussian  puff  model  for  C  (see  Ref.  19): 


C  = 


Q 


(2tt)3/2( 
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rz  u  ry  u  rz 
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where  Q  is  the  pollutant  mass  released,  U  is  the  mean  wind  speed,  t  is  the  travel  time, 
fee  is  the  effective  puff  or  cloud  height,  and  arx ,  crry ,  and  crrz  are  the  puff  standard 
deviations  in  the  x,  y ,  and  z  directions,  respectively.  Here,  he  =  hs  +  Ah  where  h,  is  the 
source  height  and  Afe  is  the  cloud  rise  due  to  buoyancy,  x  is  the  distance  in  the  mean 
wind  direction,  y  is  the  crosswind  distance,  and  z  is  the  height  above  ground.  Equation 
(1)  describes  the  concentration  field  relative  to  the  puff  centroid. 

Peak  concentration.  The  peak  concentration  is  that  in  the  elevated  buoyant  puff  which 
could  be  carried  to  the  surface  by  a  strong  downdraft  in  the  PBL,  especially  in  the  CBL. 
The  puff  centroid  concentration  Cc  is 


Cc  = 


Q 

(27r)3/2<rricrrycrrz 


(2) 


For  simplicity,  the  puff  cam  be  considered  to  be  isotropic:  arx  =  ary  =  arz  —  <jr.  For  a 
buoyant  puff,  ar  is  proportional  to  the  puff  radius  as  discussed  below. 

If  Cc  is  used  as  an  estimate  of  the  pealc  surface  concentration,  an  estimate  should  be 
given  of  its  probability  of  occurrence  there.  One  possible  way  of  doing  this  is  to  consider 
random  puff  trajectories  due  to  the  random  vertical  velocity  w  in  the  PBL: 


zp  —  fe«  H — — h  Afe(x)  , 


(3) 
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where  zp  is  the  random  puff  height.  The  probability  P(z  <  zt)  that  the  centroid  could 
be  carried  to  the  surface  is  found  from 

P{z  <zp)=  /  P(zp)dzP  >  (4) 

Jo 


where  p(zp )  is  the  probability  density  function  (p.d.f.)  of  zp  and  zt  is  a  small  height 
near  the  surface,  e.g.,  zt  ~  <xr/2.  The  p{zp)  can  be  found  from  the  p.d.f.  of  w  [p^C^O] 
according  to10 


p(zp)  =p[™(zp)>x] 


dw 

dzp 


(5) 


where  w  and  dw/dzp  are  found  from  Eq.  (3). 

Mean  concentration.  The  mean  concentration  at  a  given  height  due  to  all  of  the  random 
updrafts  and  downdrafts  is  given  by  Eq.  (1)  but  with  the  <jrz  replaced  by  <xz,  which 
corresponds  to  the  absolute  dispersion  (i.e.,  from  Taylor’s  theory,  Eq.  17  below).  This 
mean  concentration  is  relevant  for  the  SBL  or  in  the  limit  of  a  neutral  boundary  layer 
where  a  Gaussian  w  p.d.f.  is  applicable.  However,  for  the  CBL,  a  positively-skewed  w 
p.d.f.  is  more  consistent  with  laboratory  and  field  data  and  should  be  adopted. 


For  the  CBL,  a  good  approximation  to  the  w  p.d.f.  (p™)  has  been  shown  to  be  the 
superposition  of  two  Gaussian  distributions10 
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where  Ai  and  A2  axe  weighting  coefficients  for  the  distributions  with  Ai  +  A2  —  1*  The 
wj  and  crj  (j  —  1,2)  axe  the  mean  vertical  velocity  and  standard  deviation  for  each 
distribution  and  are  assumed  to  be  proportional  to  crw.  The  xv\2W2,<Ti,ct2,  Ai,  A2  are 
found  as  a  function  of  crw,  the  vertical  velocity  skewness  S  =  w3 /  cr^  where  w3  is  the 
third  moment  of  w,  and  a  parameter  R  =  <T] /w\  =  —  <72/ W2  (Ref.  20).  This  requires  cr^,, 
which  is  parameterized  in  terms  of  to*  and  the  friction  velocity  u*  (see  Ref.  10),  and  w3 
which  is  taken  as  w3  =  0.125u;2  in  the  upper  90%  of  the  CBL. 

The  vertical  concentration  distribution  is  derived  from  pw  following  the  same  approach 
as  applied  to  continuous  plumes.10  The  resulting  expression  for  C  is  given  by 


C  = 


Q 


exp 


(2tt)3/2  (Trxary  y  —  rz 

Ai  /  (z  -  he  -  ZT)2 

- exp  - 

O’zl  V 


(  (g  —  ut)2  _  _j/M 

V  2  <r2x  2o*J 


X 


2^2a 


A2 

H - exp 

&z2 


(z  —  he  —  Z2 )' 

2^22 


:) 


where 


ajx  j  — 

cr  zj  —  jj  axid  Zj 


WjX 

IT 


with  j  —  1,2. 


(7) 

(8) 
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The  A j,  <7j ,  and  wj  (j  —  1,2)  are  the  parameters  appearing  in  Eq.  (6).  Equation  (7) 
applies  only  for  small  x  such  that  the  plume  interaction  with  the  ground  or  elevated 
inversion  is  weak. 


More  complete  expressions  for  C  corresponding  to  (4)  are  applicable  in  the  case  of 
multiple  cloud  reflections  at  the  ground  and  PBL  top. 


Dosage.  There  are  practical  advantages  in  modeling  the  dosage  when  analyzing  the  air 
quality  impact  due  to  instantaneous  sources.  The  partial  dosage  ip  is  defined  by 


ip(x,y,z,t)=  [  C(x,y,z,t')dt'  (9) 

Jo 

and  the  total  dosage  by  ip^  =  ip(x,y,z,  oo).  One  advantage  is  that  ip  should  be  a  more 
stable  statistic  than  the  concentration  due  to  the  time  integration,  and  this  has  value  in 
the  analysis  of  field  data  and  model  evaluation.  Second,  the  time-averaged  concentration 
over  1  hr  periods  or  longer  is  necessary. 


For  clouds  with  short  passage  times  over  a  receptor,  the  average  concentration  C  can  be 
obtained  from 


ip{t2)  -  ip{ti) 


(10) 


where  the'averaging  time  Ta  —  t2  —  t\.  The  puff  passage  time  is  ~  AcrTX/U  and  if  this  is 
less  than  Ta ,  then  C  =  ipoo/Ta. 


The  integration  in  Eq.  (9)  can  be  carried  out  analytically  for  limiting  forms  of  the  arx , 
Cry ,  crrz,  and  Ah  variation  with  t.  For  example,  this  can  be  done  for  Ah  —  0  and 
&rx  j  &ry  y  az  oc  t  or  oc  fly/2 .  These  and  other  forms  or  combinations  of  them  must  be 
examined  to  determine  which  of  the  physically  meaningful  cases  result  in  an  analytical 
integration  in  (9).  Otherwise,  a  numerical  integration  of  Eq.  (9)  is  necessary. 

If  it  is  assumed  that  ary  and  az  are  constant  during  the  puff  passage  time  over  a 
receptor  (i.e.,  the  passage  time  is  short),  then 


^fcoo  - 


Q 

2tvU  a  ry  a  z 


exp 


(z~h>)2\ 

/ 


(ii) 


as  pointed  out  by  Gifford.19  Equation  (11)  is  of  the  same  form  as  the  expression  for  C 
due  to  a  continuous  point  source,  but  here  Q  is  the  total  contaminant  mass  and  not  the 
release  rate. 


Cloud  rise  and  inversion  penetration.  For  neutral  air,  the  governing  equations  for 
puffs  or  thermals  give  the  thermal  rise  as  a  function  of  time  and  the  initial  momentum 
and  buoyancy,  but  experiments  must  be  conducted  to  determine  an  entrainment 
coefficient  (see  Refs.  21,  22,  23).  From  a  combination  of  theory  and  laboratory 
experiments,  Scorer22  obtained  the  following  expression  for  the  rise 

Ah  =  2.35 (MTt  +  FTt2)1/4  ,  (12) 
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where  M t  and  Ft  are  the  initial  momentum  and  buoyancy  of  the  thermal.  They  are 
given  by 

Mt  =  -^-T'lwo  and  Ft  =  ^^Tf-  >  (13) 

3  °  CpPaTa 

where  wQ ,  r0,  and  Qt  axe  the  initial  velocity,  radius,  and  heat  content  of  the  thermal, 
g  is  the  gravitational  acceleration,  cp  is  the  specific  heat  of  air,  and  pa  and  Ta  are  the 
ambient  air  density  and  temperature. 

Scorer  also  reported  that  the  puff  radius  r  was  on  average  given  by  r  =  0.25A/i*,  where 
A ht  is  the  cloud  top  height.  However,  there  was  considerable  variability  in  the  above 
coefficient  which  ranged  from  0.14  to  0.5.  The  relative  dispersion  crr  oc  r. 

Using  field  observations  from  small  munitions  and  larger  detonations,  Weil24  confirmed 
that  Eq.  (12)  was  a  good  fit  to  data  over  a  wide  range  of  times  following  the  release. 
Thus,  Eq.  (12)  is  suitable  for  the  initial  rise  of  a  cloud,  i.e.,  before  it  is  limited  by  stable 
stratification.  The  initial  heat  content  Q T  of  the  cloud  can  be  determined  from  the  mass 
of  the  detonation  using  the  conversion  1100  kcal/kg  TNT  (see  Church25). 

In  stable  air,  the  maximum  cloud  rise  was  found  by  Morton  et  al.23  to  be 


Ah  =  2.66 


(14) 


.  t* 

where  N  is  the  Brunt- Vais  alia  frequency;  N2  =  (gl&){d®/dz)  where  0  is  the  ambient 
potential  temperature. 


For  thermal  or  cloud  penetration  of  an  elevated  density  jump,  results  have  been  obtained 
from  laboratory  experiments  in  a  nonturbulent  environment.  Saunders26  derived  the 
cloud  height  history  and  maximum  penetration  height  as  a  function  of  Ft  ,  the  density 
jump  A  pi,  and  its  height.  Richards27  obtained  an  empirical  expression  for  the  fraction  P 
of  the  cloud  penetrating  the  jump: 


P  ~  1 


-  0.5 


A  pi 

ApTi  ’ 


(15) 


where  A pTi  Is  the  average  density  excess  of  the* cloud  when  it  reaches  the  density  jump. 
The  A  pTi  can  be  estimated  from  Ft  and  the  cloud  radius  (r)  growth  law,  r  oc  Ah. 

Initial  criteria  for  the  cloud  fraction  1  —  P  trapped  below  a  thin  inversion  can  be 
developed  from  the  above  results.  This  can  then  be  used  in  the  dispersion  model. 
However,  the  problem  should  be  pursued  further  to:  1)  develop  consistency  between 
the  approach  used  for  clouds  or  thermals  and  those  used  for  plumes  (e.g.,  Briggs,28 
Mamins,29  and  Weil30),  2)  extend  the  model  for  P  to  a  thick  elevated  inversion 
characterized  by  the  d®  jdz,  and  3)  conduct  further  laboratory  experiments  on  thermal 
penetration  of  density  jumps  and  thick  inversions.  The  latter  experiments  should 
be  conducted  in  both  the  presence  and  absence  of  convective  turbulence  below  the 
inversion. 
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Dispersion  parameters.  Puff  or  relative  dispersion.  For  detonations  (instantaneous 
sources),  the  puff  growth  is  initially  dominated  by  buoyancy-induced  entrainment  and  r 
follows  r  oc  Ah  oc  *1/2 

as  given  above.  The  puff  should  also  grow  due  to  the  ambient 
turbulence  in  the  inertial  subrange  although  *he  observational  base  for  this  (in  the  case 


c  loyant  sources)  is  not  well  defined.  B?  on  modeling  for  plumes,8’31  a  tentative 


expression  for  the  puff  or  cloud  radius  gro  is 


dr 

—  =  «!  ,  +  a2ve  , 


(16) 


where  wp  is  the  puff  rise  velocity,  ve  =  (2 er)1/3  is  an  inertial-range  velocity,  e  is 
the  ambient  turbulent  energy  dissipation  rate,  and  <*1,0:2  are  empirical  entrainment 
coefficients.  Equation  (16)  represents  a  simple  superposition  of  the  entrainment  due  to 
buoyancy-induced  turbulence  and  ambient  inertial-range  turbulence. 

The  analogous  expression  (to  16)  for  plumes  was  used  recently  to  model  the  mean  and 
rms  fluctuating  concentrations  due  to  a  buoyant  plume  in  the  CBL.  The  approach 
produced  fair  agreement  with  the  Deardorff  and  Willis  laboratory  measurements  in  a 
convection  tank.32 


Equation  (16)  is  a  first  attempt  at  a  difficult  problem  and  one  where  laboratory  and  field 
data  would"  be  invaluable.  In  particular,  convection  tank  measurments  of  buoyant  puff 
dispersion  and  concentration  fields  would  be  very  beneficial  to  the  modeling  program. 

Absolute  dispersion.  For  a  sufficiently  long-duration  burn  (to  be  defined),  the 
“instantaneous”  or  short-time  averaged  concentration  could  be  determined  from  a 
plume  model  (Gaussian  or  p.d.f.  approach)  with  absolute  dispersion  for  az  (Eq.  17)  and 
relative  dispersion  (<rrv)  for  the  lateral  component.  However,  longer  time  averages  (say 
>  several  minutes)  can  be  determined  using  absolute  or  plume  dispersion  parameters 
for  both  the  y  and  z  components.  The  lateral  (<xy)  and  vertical  ( crz )  plume  standard 
deviations  can  be  found  from  a  parameterization  of  Taylor’s  statistical  theory:15 


crvt 

(T 


(1  +  (/27>,*)>. 


where  Tiy,Ti,z  axe  the  Lagrangian  time  scales  for  the  v  and  w  components  and  t  =  x/U. 

The  time  scales  can  be  parameterized  using  expressions  such  as  Tlz  oc  aw/h , 

Tlz  oc  <T^,/e,  etc.  as  done  previously.14’15’20  The  PBL  variables  necessary  in  these 
and  other  expressions — cru,av,aw,e,  surface  fluxes,  etc. — would  be  obtained  from  the 
mete-  logic'  platform. 

Compiex  terrain.  The  treatment  of  dispersion  in  hilly  or  complex  terrain  will  vary 
depending  on  the  wind  field  input.  In  the  case  of  wind  profile  measurements  only  (no 
diagnostic  or  prognostic  modeling),  the  dispersion  model  focus  would  be  on  the  cloud 
impaction  about  the  windward  side  of  a  hill.  The  approach  would  be  similar  to  that 
used  in  the  EPA  CTDMPLUS  model33’34  or  a  simpler  method.16  This  accounts  for  flow 
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speedup  over  a  hill,  plume  deformation,  turbulence  changes,  and  their  effects  on  the 
surface  concentration  through  a  modification  to  the  Gaussian  plume  model.  In  addition, 
it  accounts  for  the  concept  of  a  dividing  streamline  height  ( Hc )  in  stably-stratified  flow, 
where  ambient  air  tends  to  travel  around  a  hill  for  z  <  Hc  and  over  the  hill  for  z  >  Hc . 
This  approach  accounts  for  dispersion  about  the  first  hill  downwind  of  a  source  and  has 
obvious  limitations  for  sources  in  complex  terrain  consisting  of  many  ridges,  hills,  and 
valleys. 

Wind  field.  There  are  three  general  categories  of  wind  field  input  to  the  puff  model 
that  are  being  considered.  As  noted  below,  these  would  be  used  differently. 

1.  Observed  vertical  profiles  of  the  time- varying  wind  at  a  single  (x<y)  location. 
These  are  obtained  from  the  mobile  meteorological  platform.  For  modeling,  the 
observed  winds  would  be  considered  representative  of  the  wind  field  over  some 
short  range  (perhaps  10  to  20  km  or  so);  obviously,  this  depends  on  the  terrain. 
The  puff  displacements  would  be  tracked  using  the  wind  components  for  each 
sequential  time  interval.  This  would  be  used  in  the  most  routine  applications  and 
for  assessing  air  quality  impact  with  historical  meteorological  data  as  input. 

2.  Diagnostic  wind  model.  This  approach  uses  observed  winds  over  an  xyy  domain 
(mobile  platform  and  other  data)  coupled  with  the  continuity  equation  and  an 
interpolation  scheme  to  obtain  a  mass-consistent  wind  field  (e.g.,  Ref.  35).  This 

,  could  only  be  practicable  at  sites  where  adequate  wind  measurements  (a  grid)  are 
available  and  probably  only  for  selected  meteorological  scenarios;  e.g.,  this  would 
not  be  used  with  historical  wind  data  for  every  hour  of  the  year.  This  approach 
has  particular  limitations  in  complex  terrain  and  for  stably-stratified  flow. 

3.  Prognostic  wind  model.  This  approach36  could  be  used  for  selected  meteorological 
scenarios  with  observed  winds  from  the  mobile  platform  and  other  sources  as 
input.  When  used  with  four  dimensional  data  assimilation  (e.g.,  Refs.  37  and 
38),  this  could  be  the  most  general  approach  for  obtaining  the  wind  field.  Key 
limitations  are  the  computational  resources  necessary  and  the  grid  resolution. 

Turbulence  field.  The  profiles  of  cru^(Tv^aw  would  be  obtained  from  combined  use  of 
the  observed  <rw  profiles  (lowest  200  m),  the  <rU}av^aw  surface  observations,  surface  heat 
and  momentum  flux  measurements,  and  parameterizations  of  the  turbulence  variables 
at  other  sites  (similarity  profiles,  e.g.,  see  Refs.  9,  10,  15).  The  parameterizations 
would  be  used  for  guidance  as  site-specific  turbulence  profiles  may  be  generated;  a 
parameterization  of  e  also  would  be  needed.  In  addition,  the  PBL  depth  would  be 
obtained  from  the  mini-lidar. 

Other  effects.  The  applications  model  also  will  contain  descriptions  and  expressions 
for  dealing  with  dust  generated  by  the  cloud,  deposition  of  particles  (e.g.,  Ref.  39),  and 
possibly  concentration  fluctuations. 
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Research  Model 

A  Lagrangian  statistical  model  is  being  considered  for  addressing  several  aspects  of 
short-range  dispersion.  In  this  approach,20  one  follows  “particles”  in  a  turbulent  flow 
given  1)  the  Eulerian  velocity  statistics,  or  2)  the  time-dependent  Eulerian  velocity 
fields;  the  latter  are  obtained  from  large-eddy  simulations  (LES)  of  the  PBL.  Currently, 
the  Lagrangian  approach  is  being  used  to  model  the  fluctuating  as  well  as  the  mean 
concentration  field  due  to  a  passive  scalar  source  in  the  CBL,  using  LES-generated 
velocity  fields. 

The  dispersion  of  buoyant  plumes  has  been  computed  using  a  hybrid  Lagrangian  model 
that  relies  on  parameterized  profiles  of  the  Eulerian  velocity  statistics.32  This  deeds 
with  both  the  fluctuating  and  mean  concentration  distributions.  The  modeling  can  be 
extended  to:  1)  treat  buoyant  puffs,  and  2)  use  the  LES  fields  as  input  rather  than  the 
parameterized  turbulence. 


MODEL  EVALUATION 


At  present,  the  OB/OD  dispersion  model  development  includes  plans  for  testing  the 
model  with  three  types  of  data  bases. 

1)  Laboratory  data.  As  noted  earlier,  laboratory  experiments  on  cloud  or  thermal 
penetration  of  elevated  inversions  would  be  useful  for  model  testing.  We  plan  to 
conduct  such  experiments  in  a  salt-stratified  tank  in  the  absence  of  any  ambient 
or  background  turbulence.  The  experiments  would  be  similar  to  those  conducted 
by  Saunders26  and  Richards27  with  turbulent  thermals  except  that  we  will  include 
a  constant  density  stratification  above  a  neutral  layer  in  addition  to  a  density 
jump  above  the  neutral  layer.  In  addition,  it  is  planned  to  test  ambient  turbulent 
dispersion  aspects  of  the  model  using  instantaneous  releases  in  a  laboratory 
convection  tank.  The  experiments  are  planned  for  the  EPA  Fluid  Modeling 
Facility  in  Research  Triangle  Park,  NC. 

2)  Existing  field  data.  A  survey  of  existing  data  bases  on  the  rise  of  buoyant  clouds 
and  short-duration  plumes  from  surface  releases  will  be  made.  It  is  anticipated 
that  such  such  data  exist  at  military  installations  with  that  used  by  Weil24  from 
the  White  Sands  Missile  Range  as  an  example.  The  latter  included  cloud  rise 
from  detonations  as  well  as  a  plume  from  an  oil  fire.  These  data  would  be  used 
to  test  the  initial  buoyant  rise  phase  of  clouds  and  plumes  in  the  PBL. 

3)  Future  field  experiments.  Currently,  it  is  planned  to  conduct  future  field 
experiments  on  the  rise  and  dispersion  of  buoyant  clouds  and  plumes  from 
OB/OD  sources.  In  addition  to  detailed  meteorological  data  from  towers  and  the 
mobile  platform,  we  will  track  and  measure  the  dispersion  of  the  airborne  clouds 
and  plumes  with  a  lidar,  i.e.,  to  obtain  the  cloud  geometry.  We  also  will  consider 
ambient  concentration  measurements  of  cloud  constituents  with  the  feasibility  of 
such  measurements  determined  by  model  calculations  and  instrument  capabilities. 
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